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Quantum-mechanical systems having two discrete energy levels are ubiquitous in nature. For 
crossing energy levels, depending on how fast they approach each other, there is a possibility 
of a transition between them. This phenomenon is known as Landau-Zener tunneling ^"^ 
and it forms the physical basis of the Zener diode, for example. The traditional treatment 
of the Landau-Zener tunneling, however, ignores quantum-mechanical interference. Here 
we report an observation of phase-sensitive interference between consecutive Landau-Zener 
tunneling attempts in an artificial two-level system formed by a Cooper- pair- box qubit 
We interpret the experiment in terms of a multi-pass analog to the well-known optical Mach- 
Zehnder interferometer. In our case, the beam splitting occurs by Landau-Zener tunneling at 
the charge degeneracy, while the arms of the Mach-Zehnder interferometer in energy space 
are represented by the ground and excited state. Our Landau-Zener interferometer can be 
used as a high-resolution detector for phase and charge owing to interferometric sensitivity- 
enhancement. The findings also demonstrate new methods for qubit manipulations. 

Landau-Zener (LZ) tunneling is a celebrated quantum-mechanical phenomenon, taking place 
at the intersection of two energy levels that repel each other due to a weak interaction ^. The LZ 
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theory, developed in the early 1930's in the context of slow atomic collisions ^"^ and spin dynam- 
ics in time-dependent fields demonstrated that transitions are possible between two approaching 
levels as a control parameter is swept across the point of minimum energy separation. The LZ tun- 
neling is often used as a tool for determining level separations, for example, in molecular NMR ^. 
The probability of a Landau-Zener tunneling transition is given by 



where v = d{ei — So)/dt denotes the variation rate of the energy spacing for noninteracting levels, 
and 2A is the minimal energy gap. 

Yet for quantum-mechanical systems, more fundamental is the transition amplitude, which 
allows one to describe interference. As two colliding atoms approach each other, their electronic 
levels may cross. The probability amplitudes evolve along either of the two potential curves and 
may interfere, when the levels cross again after the collision. The wave-function phase accumu- 
lated between the incoming and outgoing traversals varies with the collision energy giving rise 
to Stueckelberg oscillations, observed in atomic systems in the populations. Typically, how- 
ever, the phase is large and rapidly varies with energy, which allows one to average over these fast 
oscillations neglecting the interference. 

Recently, quantum coherence in mesoscopic Josephson tunnel junctions has been investi- 
gated extensively 6.11.12^ since they might provide a realistic platform for quantum-information 
processing. In these artificial quantum systems, energy scales can easily be tuned into a range 
feasible for study of fundamental phenomena. We used a system of such mesoscopic junctions 
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to obtain the first evidence of quantum interference associated to Landau-Zener tunneling in non- 
atomic systems. 

We used a charge qubit based on a Cooper-pair box (CPB) that we turned into an analog of 
the optical Mach-Zehnder interferometer. In this device, a beam is split into two partial waves, 
which interfere after a single passage through the system. In our case, the LZ tunneling provides 
the mechanism of the beam splitting and occurs when the gate charge of the Cooper-pair box is 
swept across the degeneracy (Fig. [T]). The split beams follow the ground and excited states of 
the CPB and recombine at the subsequent degeneracy point. We find a collection of different 
types of interference patterns which can be described using the basic principles of Mach-Zehnder 
interferometers. Our interferometric observations are made possible by the non-invasive character 
of our dispersive measurement method 

Our superconducting Mach-Zehnder interferometer is made using a single-Cooper-pair tran- 
sistor (SCPT) embedded into a small superconducting loop (Fig.|2l). An SCPT consists of a meso- 
scopic island having capacitance C, two small Josephson junctions, and of a nearby gate electrode 
used to polarize the island with the (reduced) gate charge Ug = CgVg/e. The island has the charg- 
ing energy Ec = e^/ (2C) ~ 1 Kelvin, and the junctions have the generally unequal Josephson 
energies Ej(l ± d), where d quantifies the asymmetry. With d = 0, SCPT Hamiltonian is then 
Ec{n — UgY — 2Ej cos (0/2) cos(6') — CgVg/2. Here, the number h of extra electron charges on 
the island is the quantum conjugate variable to 9/2, where 9 is the superconducting phase on the 
island. The SCPT is then equivalent to a Cooper-pair box, but with an effective Josephson energy 
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of 2i?jcos(0/2) tunable by the superconducting phase across the two junctions, = 27r$/$o- 
Here, $0 = h/(2e) is the magnetic flux quantum. Diagonalization of the Hamiltonian leads to 
energy bands Ek{ng, 0) (supplementary information on bands: see Ref. '^). 



When Ec ^ Ej, the Hamiltonian is conveniently written in the eigenbasis {|n)} of the 
charge operator h, taking only two charge states into account. Then the Hamiltonian of the CPB 
becomes 



^ dun) -A ^ 
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= --B^a^ - -B^a^ (3) 
where e= \B, = 2Ec{l - Ug) and A = ^B^ = Ej cos{(f)/2). 

We display here the Hamiltonian both in the matrix form Q and in the spin form, Eq. (jSj). 
The asymmetry d 7^ in Josephson energies would limit the minimum value for the off-diagonal 
coupling I A|. The eigenvalues Eo{ng, 0) and Ei{ng, 0) are the two lowest bands as illustrated by 
Fig. [It- By |0) and we denote the corresponding wave functions. 

We analyze the level crossing process using the energy diagram in Fig.[lt. As Ug is lowered 
and then increased (similar to the interatomic distance during a collision), the system can follow 
either of two possible paths: AOCOD and AOBOD. The probability to follow either path is a 
product of two independent events: Plz for making a transition and 1 — Plz for staying on the 
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same level, which gives the total transition probability 

Pad = 2Plz(1 - Plz) • 
In this treatment interference between tunnel attempts has been neglected. 
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Subsequent LZ tunneling events with time interval Tp can interfere, provided phase coherence 
is preserved and these events do not overlap < Tp < Tcoh- Here, the time of an LZ-tunneling 

event is ~ ■ max(l, /S? /hv). In charge qubits, it is easy to make ^ Tcoh where 

the coherence time is Tcoh = niin(Ti,T2) with Ti and T2 corresponding to the relaxation and 
dephasing time, respectively. For example, by taking A = 2 GHz and f = 40 GHz per 1 ns, 
we obtain ~ 0.1 ns, which is well within experimental reach. The interference of consecutive 
tunneling attempts can be viewed as two partial waves, describing the propagation along either the 
lowest band or the first excited band. This is similar to an optical Mach-Zehnder interferometer as 
illustrated in Fig. [TJ). 



Away from the crossing region, the eigenstates |0) and |1) accumulate the dynamical phase 

1 



^ = _ y,(0) = _ I [E,{ng{t)) - Eoin.m dt . 
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In addition, each pass of the level crossing results in an LZ event with probability amplitudes given 
by (cf.Fig.Uh) 



1 , ^ ^ 
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VI - PLzexp(i(^5) 



LZ 
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i\/PvL VI - PLzexp(-z0s) 

Here, 4>s = 4>s — 7r/2, where the Stokes phase 05 depends on the adiabaticity parameter /hv 

(cf. Eq. dU)). In the adiabatic limit, (ps 0, but in the sudden limit, 05 = ti/A. 
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Let us assume a fast gate charge sweep of the form ng(t) = UgQ + Srij-^ sm{ujj.it), where the 
constant level UgQ means that the sweep is in general offset from the crossing point. One cycle 
of continuous driving in our CPB takes the system twice through the crossing point and involves 
two dynamical phase shifts cpi and cpR, on the left and right sides. For a single cycle, we add the 
amplitudes along the two branches in Fig.[lj5, to find the probability of reaching the point D: 



P, 



AD 



2Plz(1 - Plz) 1 + cos(¥Pl - 205) 



(7) 



Clearly, the maximum transition probability is reached when the total phase v^l — 2(f)s is a multiple 
of 271. Under continuous driving, one obtains a multi-pass Mach-Zehnder model. In this case it 
can be shown that the maximum population of 1 1) (constructive interference) is reached when both 
dynamical phases satisfy the condition mentioned above. 



V^L,R — 205 are multiples of 2n 



(8) 



For example, in the adiabatic limit, v^l.r have to be odd multiples of vr. The resonance conditions 
in Eq. ^ are seen overlayed in Figs. |3] and HI (see below) as the black solid and dashed lines. 

Our experimental scheme is illustrated in Fig. |2| The measurement signal tracks the time 



average, under a strong LZ drive, of the Josephson capacitance Ccfr oc 



probed at /„ 



803 MHz (see Methods). In the first approximation, the energy E here can be taken as the average 
energy stored in the qubit: E = PqEq + piEi, where the band energies are weighted by their 
average populations po and pi, respectively. Thus, we have 



Ceff = PoCls + PlCl 



cS • 



(9) 



In this way, however, we neglect all the relaxation phenomena that take place on time scales 1 / fp-. 
If the relaxation rate (Ti)~^ ^ fp then changes in pq = Po{Vg) have to be taken into account in the 
response. As will be seen below this is in fact crucial in understanding our experimental results. 

We have made extensive scans of the reflection coefficient of a CPB by varying the LZ drive 
frequency = 0.1 — 20 GHz, and its amplitude Sn^f = — 3 electrons, as well as the bias UgQ 
and (p. The Josephson capacitance deduced from the phase shift of the reflected wave arg(r) at 
/rf = 4 GHz when changing Suj-f and n^o is illustrated in Fig. 13 We observe a clear interference 
pattern whose main characteristics agree with those expected for coherent LZ tunneling: 1) There 
is an onset of the interference speckles with a distinct value for Ugo where the rf drive just reaches 
the avoided crossing, linearly dependent on the AC drive amplitude; 2) The density of the dots is 
proportional to 1/ /rf in the direction of Ugo as well as 5nj-{, 3) The central part of the interference 
patterns displays the curved dot rows, in a similar fashion as in the overlayed patterns, 4) the pattern 
loses its contrast at a certain value towards lowering drive frequency, here at /rf ~ 2 GHz. Note 
also that there are destructive interference dots at high drives, where the qubit remains basically on 
the lowest level (cf. the "coherent destruction of tunneling" ^'^). 

The periodicity in UgQ is clearly 2e at low levels of rf-excitation. At excitations on the order 
of e/2, there is an appearance of a shifted, additional pattern, which makes the signal almost 
e-periodic. The origin of these odd sectors can be understood by looking at the energy levels 
displayed in Fig. [H When the rf-drive brings the system past a crossing point of Ei and Eq, it 
becomes energetically favorable to enter an odd particle-number state, resulting in a shift by e in 
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the interference pattern. The odd states appear to be rather stable , and the contrast of the "odd 
sectors" is almost as strong (see Fig.|31l. 



According to Eq. (HJ, the phase difference, 



4Ec{l-ngo) 



(10) 



is a multiple of 2n at resonances, implying the location of the population peaks on the lines of 
fixed gate bias Ugo with spacings Augo = hwrf/ {2Ec) (we approximated the dynamical phase by 
that for non-interacting levels). This linear dependence of the spacings on frequency is illustrated 
in Fig. 13;. The fitted line yields Ec = 1.1 K, which is about 25 % higher than we obtained from 
the rf-spectroscopy 

As usual, interference effects are prone to decoherence and our interferograms are suitable 
for studying dephasing and relaxation in qubits as proposed by Shytov et al. ^. The suppression 
of interference in our data at low /rf is due to the loss of phase memory over a single LZ cycle. 
Indeed, phase fluctuations suppress the contrast of oscillations in Eq. © (see Ref. for a more 
detailed analysis). This way, we find a quick estimate, averaged over Ug, for the coherence time of 
our qubit Tcoh ~ 0.5 ns. 

To account for decoherence in a detailed manner, we solved the phenomenological Bloch 
equations ^^'^^ which describe the dynamics of the magnetization M = (S) of a pseudospin-1/2 
(a two-level system): 




B X M 



1 



1 



(11) 
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Here, the pseudomagnetic field is given by Eq. ©. The parameters Ti and T2 describe the relax- 
ation of the 2;-component of magnetization towards equilibrium and the relaxation of transverse 
magnetization to zero, respectively. Assuming that decoherence is dominated by charge noise, we 
write 



Here, the angle tan r] = B^/ describes the dependence on the gate bias, and the voltage fluctua- 
tions are involved via S'x(co') = (^26^) S'v'(co'). For Ohmic dissipation, 5'x(t^) = 27r/l^at<; coth ; 



where the coefficient a characterizes both the bath and its coupling to the qubit: a = (1 + [Ci + 
C2\/ Cg)^'^'^R. For our sample, a ~ 10^^ due to strong coupling to the environment via the high 
gate capacitance and parasitic capacitance in the resonator inductance. 

The results of the simulations for a strong driving Bz = 2Ec{l — rigo + Suj-f sin(u;i.ft)) are 
illustrated in Fig. HI It displays the population of the ground state on the UgQ vs. 5nj.{ plane for 
Ohmic bath with a = 0.04. The variation of the level population agrees well to the overlay ed inter- 
ference pattern obtained using the Mach-Zehnder model, except for the regions near the edges (the 
inclined white lines in Fig.©. The additional structure in this region is a signature of multiphoton 
transitions: Under the condition of the n-photon resonance (4£'c(l — ^^^) = ^ ' ^rf) the effective 
coupling in the pseudospin rotating frame is A = A • Jn{'^Ec5n^f / fuJrf), and from the stationary 
solution of the Bloch equations (with gate-independent relaxation terms) we find the population 
response oc TiT2A'^/{l + TiTgA^ + {\E^{ngo - n^^^)T2f). The Bessel functions describe the 
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onset of the resonance as well as the interference pattern at higher drives. This expression indicates 
that destructive interference occurs at amplitudes corresponding to zeros of the Bessel functions, 
and that the destructive interference spots are sharper, in accordance with the data. 

The range of the capacitance variation on the lower level, cj,^ ~ 0.2-1.5 fF, deduced from 
the data measured without the LZ drive, are in accordance with the ground band curvature, shown 
in Fig. 12]:. As the maximum population on the upper level can only reach 50 %, the weighted 
average, Eq. ©, should always be larger than 0.2 fp (see Fig. |21). Clearly, there are regions in 
Fig. |3l where the response looks stronger. Hence, instead of Eq. we have to calculate a time 
average (— |^), whose magnitude can be substantially increased by relaxation phenomena. 

Assuming very fast relaxation, i.e., with the populations tracking the instantaneous AC gate 
charge Ug = UgQ + Sn,.f sin(ci;i.ft), we find a capacitance of good uniformity, in a similar manner 
to that of the measured interference patterns. However, at reasonable values of parameter a, the 
swing of the capacitance is too large indicating that the assumption of instantaneous relaxation is 
too strong. 

With a finite relaxation Ti, we have calculated the spin dynamics having a weak measurement 
signal of amplitude Sn^c = CgSV^c on. Then, we can use the linear-response theory to extract the 
capacitance. We calculate the time-dependent expectation value for the effective charge (Q) (t) = 
Tr(p * Qcs), where Qes = —dE/dVg, and the density matrix is expressed in the energy eigenbasis. 
From (Q) (t) we pick up the quadrature components, Q^^.^^ and Q^ovlx ' the measurement frequency. 
Using the definition of impedance, we may solve for the capacitance C^s = .,7°"* • The 

^•^in 

10 



resulting capacitance at /rf = 4 GHz is illustrated in Fig. 13 The full swing of the capacitance 
over the UgQ vs. 5n^f plane is very sensitive to the parameter a. The values a = 0.04 and Suac = 
0.03 were taken for the calculation in order to match the measured pattern. This corresponds 
at the degeneracy point to T2 ~ 0.5 ns which is close to the estimates obtained both from the 
oscillation contrast, and microwave spectroscopy. The calculation reproduces the major features 
of the measured interferograms: 1) The spacing of the dots, 2) the size of the capacitance swing, 
3) the global uniformity of the pattern, and 4) minima and maxima on separate arcs (rather than on 
the same arc). 

Similarly as optical interferometry played a central role in the development of "photon" 
physics solid-state quantum interferometry may find many applications. We propose to apply 
the LZ interferometry for sensitive detection of phase and charge utilizing mesoscopic Josephson 
circuits 14,27,28^ where it brings about a significant increase in sensitivity. Indeed, the LZ inter- 
ferometer can be viewed as integrating phase amplifier for the superconductor phase across the 
device. The interferometer transforms tiny changes of (or magnetic flux $) into a huge modula- 
tion of the wave-function phase by basically integrating the hatched area in Fig.[l]^^, but at the 
expense of reduced measurement strength. Eventually, it amounts at least to ~ 27r-fold increase in 
detector sensitivity, with significantly reduced disturbance due to the measurement signal. Another 
possibility would be to separate the LZ modulation drive to the other quadrature, namely 0, and 
measure charge. The modulation in the off-diagonal components would result in slightly different 
interference phenomena, as predicted recently by Hanggi and coworkers 
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The coherent Landau-Zener process observed in this work also serves as a suitable method, 
intermediate between the two existing schemes, for manipulations of superconducting qubits. Thus 
far, qubit operations have been carried out either using lengthy rf-pulses which induce Rabi oscil- 
lations, or by sub-nanosecond "hard" pulses ^ shaped rectangularly in time. The LZ manipulation 
would offer an ultra-short clock period <^ ns similarly to hard pulses, but would be more precise 
due to basically single-frequency drive, and also technologically simpler. 
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Methods 

Our investigations of the LZ tunneling are based on measuring the quantum (or "Josephson") 
capacitance of a CPB ^^'^^ This capacitance is related to the curvature of band k similar to 
the effective mass of an electron in a crystal: 



The difference in the Josephson capacitance for A; = 0, 1 allows us to determine the state of the 




d^Ek{(j),ng) 



Cld^Eu{<P,ng) 



(14) 



CPB. 
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We perform low-dissipation microwave reflection measurements 15,27,28 ^ series LC res- 
onator in which the box effective capacitance, Eq. (fT?b . is a part of the total capacitance Cs + . 
The resonator is formed by a surface mount inductor of L = 160 nH. With a stray capacitance 
of C5 = 250 fF due to the fairly big lumped resonator, the resonant frequency is /o = 800 MHz 
and the quality factor is Q — 16 limited by the external Zq = 50 Q. When C^'^ varies, the phase 
arg(r) of the reflected signal Kut = ^Vin changes, which is measured by the reflection coefficient 
T = {Z — Zq)/{Z + Zq) = roe*'^'^^^'"^ Here, Z is the resonator impedance as marked in Fig.|2l 
Since we work rather far from matching conditions, the magnitude of the reflection coefficient Fq 
remains always close to one. The variation in arg(r) due to modulation in C^'^ is up to 40° in our 
measurements, corresponding to a shift of resonance frequency A/p ~ 6 MHz. In addition to band 
pass filtering, we used two circulators at 20 mK to prevent the back-action noise of our cryogenic 
low-noise amplifier from reaching the qubit. In all the measurements, the probing signal Vin was 
continuously applied. 
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Figure 1 : a, schematic view of interference of successive Landau-Zener (LZ) tunneling attempts in 
Cooper-pair box (CPB). From the initial state at A, the state evolves through the first LZ event at the 
avoided level crossing at Ug = 1 (O) either towards B (no LZ tunneling) or C (with LZ tunneling). 
After the turning points B and C, the final state D is reached either by a second LZ tunneling or 
by remaining on the excited band, respectively. The hatched area determines the dynamical phase 
shifts (pL,R- The level spacing Ei—Eq varies roughly as AEc{^ — ng), i.e., linearly as in the generic 
LZ-tunneling problem with linearly crossing energy levels. 2Ao ~ 3 GHz denotes the minimum 
gap in our experiments. The dashed line represents the lowest energy of odd parity state E'^'^'^. 
Even states with energy larger than the odd state value will always try to relax to the odd state, 
b, interpretation of the LZ interference in a as a multi-pass analog of the optical Mach-Zehnder 
interferometer. LZ events correspond to the beam splitters which divide the wave function into 
two partial waves, with the (probability) amplitudes as marked. The mirrors play the role of the 
dynamical phases (/?(°'^^ = h^^ J E^^idt picked up away from the avoided crossing. 
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Figure 2: a, schematics of our experiment. The resonant frequency /o of the lumped-element 
LC circuit is tuned by the Josephson capacitance Ccs of the CPB shown in the SEM image. The 
junctions of the split Cooper-pair box had an average resistance of i? = 23 kTl each, corresponding 
a maximum Josephson energy of the box 2Ej = 12.5 GHz, which could be tuned down to 2.7 
GHz by magnetic flux The capacitance the junctions amounts to Ci + C2 ~ 0.34 fF, yielding a 
Coulomb energy of e^/2(Ci + C2 + Cg) = 1.1 K. b, illustration of the phase shift arg(r) of the 
reflected microwaves at a fixed measurement frequency fm while Ccs increases. C, the Josephson 
capacitance calculated for the two lowest levels of our CPB with Ej/Ec = 0.27 and asymmetry 
d = 0.22 at = 0. Note that variations from the "classical" capacitance level, (l/Cg + l/(Ci + 
6*2))^^ — 0.2 fp, are opposite for the two levels: ground level capacitance C^g > while the 
excited state has Cgg < 0. 
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Figure 3: a, interference pattern obtained from the measured microwave phase shift at = 
4 GHz. The phase bias = 0, corresponding to the level repulsion of 2A = 2Ej = 12.5 
GHz. The color codes indicate the equivalent capacitance obtained using standard circuit formulas 
(see Methods). Around Ugo = —1, the imposed lines illustrate the conditions of constructive 
Landau-Zener interference, which is expected when the phases picked up to the left of = — 1, 
— vr + 205 (solid), and to the right, — tt + 205 (dashed), are integer multiples of 2n (see 
Eq. (HI). Here, due to the almost adiabatic limit, the Stokes phase 05 = 0. The population of the 
upper state is expected to be the strongest (red) when both conditions are satisfied. Around Ug = 1, 
is imposed the equicapacitance contour for Ccs = fF obtained from the spin dynamics simulation 
using Bloch equations (compare to Fig.|5l) which agrees quite well with both the resonance grid as 
well as with the data. The interferogram b is similar to a but with /rf = 7 GHz, and 05 = 7r/4 due 
to operation in almost the "sudden" limit of LZ tunneling. C, the averaged spacing of the central 
interference peaks in gate offset (as depicted in a), with the phase bias values of (squares) and tt 
(circles). The expected linear behavior yields a fit Ec = 1.1 K. 
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Figure 4: Interference-induced variation in level populations (color codes indicate relative popu- 
lation Po of the ground state) obtained from the simulation of Bloch equations for the qubit using 
a driving field of Bz{t) = 2Ec{l — %o + ^^rf sin(co'rft)), with uj^{/{2n) = 4 GHz. Ohmic bath 
was assumed for the dissipation and the parameter a of Eq. (fT^ was set to 0.04. The inclined 
white lines illustrate the locations of reaching the degeneracy point during the sweep for the first 
time when ramping up 5nj-{ at fixed Ugo (starting point with Suj-f < 1 — Ugo). The solid and dashed 
black lines indicate the locations where conditions for 27r-multiple phase shifts on the both sides 
of the degeneracy point are fulfilled, Eq. (ISJ), with $5 = (similarly as in the experimental data in 
Fig. 13). 
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Figure 5: Josephson capacitance Ccs, calculated as arising from curvatures of the two lowest bands 
and interband relaxation, using Bloch equations and linear response theory. The dissipation param- 
eter a of Eq. (fT2t was set to 0.04 and the amplitude of the ac-excitation at 803 MHz corresponds to 
0.06e peak-to-peak. The capacitance variation is intermediate between simple averaging and the 
fast relaxation approaches, and it agrees quite well with the measured results in Fig. 12 the compar- 
ison with data is performed by the displayed contour graph which describes the equicapacitance 
curve for Ccs = fF obtained from the simulation. 
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